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Abstract. In the spirit of the White-Bear version of fundamental measure theory 
we derive a new density functional for hard-sphere mixtures which is based on a 
recent mixture extension of the Garnahan-Starling equation of state. In addition 
to the capability to predict inhomogeneous density distributions very accurately, 
like the original White-Bear version, the new functional improves upon consistency 
with an exact scaled-particle theory relation in the case of the pure fluid. We 
examine consistency in detail within the context of morphological thermodynamics. 
Interestingly, for the pure fluid the degree of consistency of the new version is not only 
higher than for the original White-Bear version but also higher than for Rosenfeld's 
original fundamental measure theory. 



1. Introduction 

In the mid 1970s, density functional theory, which was originally formulated for quantum 
systems, has been extended to systems that follow classical statistical mechanics p. 
Since then, density functional theory of classical systems (DFT) has developed to an 
indispensable tool for the study of inhomogeneous systems such as crystals, fluids in 
confined geometries [2j, liquid- vapor interfaces and wetting and drying on substrates 
(for a recent review see DFT is based on the fact that there exists a functional Q[p] 
of the spatially varying particle number density p(r) which possesses two properties: 
(i) it is minimized by the equilibrium density Po(r), and (ii) the minimum value fi[po] 
equals the grand potential fl of the system. These properties give rise to the variational 
principle 6Q[p]/6p = for p(r) = po{r). 
One can decompose Q[p] as 

m = m + J drp(r)(V;xt(r) - p) , (1) 

where p is the chemical potential, Vcxt(r) the external potential acting on the particles 
and J-'[p] is a unique functional corresponding to the intrinsic Helmholtz free energy of 
the system in equilibrium. 
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In principle, Eq. with the variational principle constitutes an excellent tool for 
(numerical) calculations of Po(r) and hence the grand potential r2[po(r)] in arbitrary 
external potentials. Unfortunately, for many systems of interest, only more or less crude 
approximations for jF[p] are known. The expression for the ideal gas, however, is known 
exactly 

= J drp(r) (ln(AV(r)) - l) , (2) 

where A is the thermal wavelength of the particles and (3 = \/{k-QT) with Boltzmann's 
constant k-Q and the temperature T. The interactions among particles are described by 
the excess (over ideal gas) free energy J^cJyp] = ^[p] — ^id[p\ which for our purposes 
can be expressed as [3J-'cx = /dr$(r), where the excess free energy density $(r) is a 
functional of p(r) . 

In this work we direct our attention to mixtures of hard-spheres. For these systems, 
the exact expression for is unknown but a number of approximations can be found in 
the literature j3]. The interest in the hard-sphere system is manifold. The hard-sphere 
system serves as a reference system for fluids with short-ranged repulsion and additional 
attractive interactions among particles. The attractive part of the potential is usually 
treated perturbatively jH]. Furthermore, colloidal suspensions with quasi hard-sphere 
interactions can be realized experimentally (see, e.g. jH]) which provides a test ground 
for predictions from the field of the purely entropically driven hard-sphere systems, like 
entropic forces [Zj, asymptotic decay of correlation functions in mixtures jHUH] etc.. 

A very successful class of excess free energy functionals is formulated within the 
framework of fundamental measure theory (FMT) introduced by Rosenfeld jlH] . While 
the original FMT has the Percus-Yevick (PY) equation of state as an output, later, 
with the setup of the White-Bear version of FMT ^3 E] , the more accurate Carnahan- 
Starling (CS) equation of state was incorporated into FMT. The resulting gain in 
precision in the structure of inhomogeneous density distribution |^ |^ and in 
thermodynamics, however, has to be paid for with a slight inconsistency appearing 
on the level of the pressure fT]: the pressure in the hard-sphere fluid obtained from 
a scaled-particle theory equation differs slightly from the underlying bulk equation of 
state. The aim of this work is to build upon the White-Bear version of FMT, using a 
new mixture formulation of the CS equation of state jH] , such that this inconsistency 
is resolved. 

The paper is organized as follows. In Sec. |21 we review Rosenfeld's derivation of 
FMT. Section ini is dedicated to the presentation of the White-Bear version of FMT and 
the derivation of the new version of FMT. In Sec. ^ we give a brief introduction to 
morphological thermodynamics of hard-sphere fluids and compare the performance of 
Rosenfeld's FMT and the original White-Bear version with that of the new functional. 
Section El contains our conclusion. 
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2. Fundamental measure theory 

In a seminal paper in 1989, Rosenfeld set up FMT, allowing him to derive his successful 
free energy functional for the hard-sphere mixture ^U]- We sketch his approach briefly 
in the following. 

We consider a //-component hard-sphere mixture with spatially varying particle 
number densities Pj(r), i = 1, . . . , i/. From the theory of diagrammatic expansions j3] 
the excess free energy functional in the low-density limit is known exactly. One finds 
that the Mayer-/ function, related to the pair interaction potential Vij{r) between the 
particles of species i and j by fij{r) = exp[—pVij{r)] — 1, plays a central role. In the case 
of hard spheres, Vij is either infinite if spheres overlap and zero otherwise. As a result 
the Mayer-/ function obtains a purely geometrical meaning: fij{r) = — G(i?j -|- Rj — r) 
where Ri and Rj are the radii of the respective species and is the Heaviside function. 

Inspired by the exact excess free energy functional of the one-dimensional hard-rod 
mixture [T5tll6j. the key-idea of FMT is the deconvolution of the Mayer-/ function fij{r) 
into a sum of products with factors depending only on one of Ri and Rj. Rosenfeld's 
deconvolution reads 

- fij{\ri - Tjl) = Uq(E) ul + ul + ui + ul- ujI^ u:^- u)2<^ u^l (3) 
with four scalar and two vectorial (weight) functions 

^(r) = e(i?,-|r|), ^(r) = 5(/2,-|r|), ^(r) = ^ - |r|) , (4) 

and ^(r) = ^(r)/(47ri?i), cj*(r) = u^ir) / {4tt R^) , and uj{{r) = uji{r)/ {An Ri). The 
convolution product ® in Eq. Q is defined by 

^1®^^ = y"dra;^(r-r,) -^^(r-rj), (5) 

where the dot ■ stands for the usual product in the case of the scalar weight functions 
and for the scalar product in the case of the vectorial weight functions. Using the weight 
functions, one can define weighted densities 

^"(^) = 5Z / dr'p^(rVa(r " r') • (6) 
i=i ^ 

The deconvolution, Eq. (jS)), can be used to express the exact low-density limit of 
the excess free energy functional: 

lim = -\Y. I drdr'p,,(r)p,(r')/.,(|r - r'|) 

h3=l (7) 

= j dr(no(r)n3(r) + ni(r)n2(r) - ni(r) ■ n2(r)) . 

This result, together with the structure of the exact excess free energy functional of 
one-dimensional hard-rod mixtures, leads to the assumption that the excess free energy 
density <I>(r) can be approximated as a function of the six weighted densities only. This 
assumption guarantees to recover the exact low density limit of <I>(r). The expression for 
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the free energy density is obtained by an extrapolation of the known low-density result 
for $(r) to higher densities using thermodynamic arguments. 

We consider the case of a homogeneous hard-sphere mixture, i.e. the density 
distributions Pi(r) = = Ni/V are constant. Ni is the number or spheres of species i 
in the volume V. The excess pressure p'^^ of a fluid mixture can be obtained from the 
excess free energy density $ via 

1=1 a 

Note that the vectorial weighted densities, which are formally included in the above 
sum, actually vanish in the uniform fluid. The ideal gas contribution to the pressure 
is jSp^"^ = J2iPiy which in terms of the weighted densities reduces to = no. Hence, 
according to thermodynamics (TD), the total pressure within FMT can be written as 

PPtd = no - $ + 2^ a^"" ■ 

a ^ 

On the other hand, there is an exact relation from scaled-particle (SP) theory (see 
PH] and references therein) between the chemical potential /Xj of species i for a very 
large sphere and the reversible work required for the creation of a cavity that can hold 
the large sphere of species i: in the limit Ri ^ oo one obtains fii/Vi — > psp, where psp 
is the total pressure of the fluid mixture and Vi = (4/3)7ri?,?. In our context, this is 
equivalent to (cf. [HHITT]) 

PPsp = t-. (10) 

Obviously, we obtain a differential equation for $ by equating the expressions for 
Ptd and psp- This equation was solved by Rosenfeld who used the ansatz 

^ = fiMriQ + f2{n3)nin2 + /3(?^3)ni ■ ns + 74(^3)^2 + /5('^3)'^2n2 ■ , (11) 

with /i, . . . ,/5 being functions of the dimensionless weighted density n^. The ansatz 
Eq. (jlip combines all multiplicative combinations of the weighted densities which 
share the dimension of $, i.e. (length)"'^. There is a solution $rf of the SPT 
differential equation. The integration constants can be fixed by the following additional 
requirements: (i) in the low-density limit, Eq. ((Zj) is recovered, (ii) for the one-component 
uniform hard-sphere fluid the correct third virial coefficient is reproduced and (iii) the 
pair direct correlation function c''^^(r) is regular for r ^ 0, which enforces the prefactor 
for the term oc 112112 ■ 112. The result is 

. . nin2-ni-n2 n| - 3n2n2 ■ n2 
*RF = -^0 In 1 - ng + + . 12 

Rosenfeld's excess free energy gives a good account for many aspects of nonuniform 
hard-sphere fluids, pure [T71[TH] or mixtures [12]. However, it does not predict freezing, 
which is actually observed for the pure hard-sphere fluid at a packing fraction rj ~ 0.494. 
This deficiency can be resolved empirically by modifying the third term in Eq. ()12p 
|2lH I2T] . or more systematically by the recipe of Tarazona j22] who introduced an 
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additional tensorial weighted density in the last term of $rf- There remains another 
factor which limits the accuracy of $rf, namely the equation of state obtained from 
Eq. or, equivalently, Eq. (fTO|l . One finds the pressure 

n no nin2 nl 

l-ng (l-ng)^ 127r(l-n3)'^ 
which is the compressibility expression from the solution of the PY integral equation 
[23] . The PY pressure is in good agreement with simulations for the pure hard-sphere 
fluid at low packing fractions but close to the freezing transition it overestimates the 
pressure by about 7%. This problem can be solved by incorporating more accurate 
equations of state within the context of Rosenfeld's FMT as an extrapolation from low 
to high densities ^3 E] • In the next section we introduce a new contribution along 
these lines. 

3. The new functional 

An empirical correction of the high-density behavior of the PY compressibility result has 
been given by Carnahan and Starling [211. The CS equation of state has subsequently 
been generalized to hard-sphere mixtures resulting in the Boublfk-Mansoori-Carnahan- 
Starling-Leland (BMCSL) pressure Pbmcsl [211121] • This equation of state can be written 
in terms of the weighted densities of a homogeneous hard-sphere mixture: 

l-ns (l-ns)^ 127r(l - Tig)'^ 

Using this fact, Pbmcsl has been incorporated into FMT ^TJ^j. This was achieved 
by solving the differential equation for the excess free energy $ which is obtained by 
equating pbmcsl and the thermodynamic expression ptd as given in Eq. Q. Note that 
for the implementation of this approach a bulk fluid mixture is considered for which 
the vectorial weighted densities vanish. Hence, the solution of the resulting differential 
equation is obtained by using the dimensional ansatz Eq. without the vectorial 
contributions. For this ansatz there is a unique solution if two additional requirements 
are made: (i) the result for $ is compatible with the low-density limit Eq. ([7|) and (ii) for 
the pure hard-sphere fluid the third virial coefficient is recovered. Unlike in Rosenfeld's 
derivation of $rf the vectorial contributions have to be incorporated at a later stage. 
In analogy to $rf the substitutions nin2 nin2 — ni ■ n2 and ^ — 3n2n2 ■ n2 
are made in $. The resulting functional $wb is called the White-Bear version of FMT 

In virtue of these substitutions, $wb has the correct low-density limit Eq. ([7|) and 
the regularity of the pair direct correlation function for r — > is guaranteed. 

The White-Bear version of FMT has been shown to inherit all the good properties 
of Rosenfeld's FMT for the description of the hard-sphere fluid and improves the 
predictions of thermodynamic quantities due to the more accurate underlying equation 
of state. This becomes particularly apparent in the contact densities at a hard wall 
which are related to the pressure via a sum rule (for a comparison with simulation data 
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see, e.g. jlll)- Furthermore, one finds that the prediction of the freezing transition of 
the pure hard-sphere system agrees very well with simulations A drawback of 

$WB is, however, that the scaled particle relation Eq. (fTO|l is violated, i.e. one finds 
that ^^^^B/^n3 7^ Pbmcsl- This is of course not surprising as the equahty ptd = Psp 
unambiguously leads to Rosenfeld's $rf, if we assume that the free energy density 
is a function of the weighted densities no, . . . ,^,3 and iii, and 112 alone. Despite this 
inconsistency of the White-Bear version, the quality of the resulting density distributions 
is high ^Hj- However, analytical results obtained from the free energy density $wb 
within the context of morphological thermodynamics (j2H] and references therein) are 
affected. 

We conclude, that there is some room for improvement with respect to the self- 
consistency of the free energy density. The basis for this improvement is a new 
generalization of the CS pressure to mixtures of hard spheres which was recently 
suggested by the authors [Hj. In terms of the weighted densities, the new equation 
of state reads 

no n,n2 (l + ^2 (l " §^3 + j^l) . . 

PPcsm = \ -r, ^ \ 777-71 V? • ^^^> 

The index CSIII refers to a hierarchy of extensions introduced in Ref. ^3], where we 
showed for binary and ternary hard-sphere mixtures that pcsiii improves upon Pbmcsl 
compared to computer simulations. Even more interestingly for the present context, 
Pcsiii was constructed such that it is consistent with the scaled-particle relation Eq. (fTO|) 
in the case of the one-component hard-sphere fluid. Note, that consistency for the 
general hard-sphere mixture within the framework of FMT would always result in the 
less accurate pressure ppy. 

By following the recipe for the derivation of the original White-Bear version pT| ll2j. 
described above, we calculate a new functional based on the pressure pcsiii 

nin2 - ni ■ n2 



*WBii = - no ln(l - ng) + (l + |n^02('^3)) ■ 



1 - ns 



with 

02(^3) = (6^3 - Ml + 6(1 - ng) ln(l - n^)) /nl = l + + 0{nl) , 

03(^3) = (6r23 - + + 6(1 - 7.3)2 ln(l - 7.3)) /(4n^) = 1 - 1^3 + O(n^) . 

The new functional is an improvement of the White Bear version of FMT, as we shall 
show in Sec. El The index WBII is chosen to indicate that the new functional is Mark 
II of the White Bear functional. 

For comparison we mention that in the above notation the original White-Bear 
functional $wb is recovered with 0^^(723) = and 

0rK) = {^l - 6713 - 6(1 - nsf HI - ns)) /{4nl) = I + ^n, + 0{nl) . (18) 

We have compared predictions of our new version of FMT with corresponding 
results obtained by the original White-Bear version for a pure hard-sphere fluid and a 
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binary mixture close to a planar hard wall. We have found that the density distributions 
resulting from numerical minimization of the functional Eq. with $wb or $wBii! 
respectively, differ very little. For the pure hard-sphere fluid, this can be expected from 
the fact that the underlying bulk equation of state is the same for both versions of 
FMT and hence the contact densities at the wall have to be identical. Comparison 
with density distributions from Monte-Carlo simulations revealed that the very small 
deviations of the DFT results from the simulation data are clearly more significant than 
the mutual deviations between the two FMT versions. We conclude that the limitations 
of FMT-based density functionals cannot be considerable pushed forward by increasing 
the quality of the underlying bulk equation of state but are rather determined by the 
structure of FMT itself, i.e. the set of weight functions which are employed and hence 
the restriction to one-center convolutions. For a discussion of this topic see Ref. [29] . 
A slight improvement from the WBII version is indeed found for the description of the 
pair direct correlation function as can be inferred from comparison with simulation data 
(not shown). 

We find the main benefit of the new functional $wbii in the context of morphological 
thermodynamics. Here, the self-consistency of $wbii on the level of the pressure, i.e. 
the equality of Ptd from Eq. Q and psp from Eq. (fTUI) in the case of the pure fluid, 
is crucial for the accuracy of analytical expressions obtained within the morphological 
theory. In the next section we give a brief introduction to the theory and show examples 
which illustrate the gain from the new functional $wbii- 

4. Morphological thermodynamics 

The morphometric approach to the grand potential of a fluid around a complexly 
shaped object B (or a fluid inside a complexly shaped container) was inspired by the 
Hadwiger theorem [20] from integral geometry jSI]. The theorem states that every 
motion-invariant, continuous and additive functional of the complexly shaped object B 
depends on the shape of B via only four geometric measures: the volume V , the surface 
area A, the integrated (over the surface area) mean and Gaussian curvature C and X, 
respectively. The latter are given as 



Here Ri and i?2 are the local principal radii of curvature on the surface of B. Note that 
X is proportional to the Euler characteristic. 

While the Hadwiger theorem is rigorous, the connection to physics is not obvious 
and cannot be proven rigorously. However, there is strong numerical evidence 
[321 I2H1 I2Z1 1221 that the solvation free energy of a convex body immersed in a 
solvent away from the critical point and away from wetting or drying transitions takes 
the form 




(19) 



/S,n = pV + aA + kC + RX . 



(20) 
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where the conjugated quantities to the geometric measures of B are thermodynamic 
coefficients depending only on the temperature, the chemical potentials and the given 
interaction between B and the fluid, and among fluid particles, but not on the (complex) 
geometry of B. The thermodynamic coefficients are p, the pressure, a the planar wall 
surface tension, and k and R two bending rigidities. 

Morphometry is obviously a very useful tool for the calculation of thermodynamic 
quantities in complex geometries as it allows one to calculate the shape-independent 
thermodynamic coefficients in simple geometry. The treatment of the actual complex 
geometry (or a set of different geometries) then only requires a straightforward 
calculation of the geometric measures. For instance, the morphometric approach has 
been applied successfully to the calculation of solvation free energies of a protein in 
various geometrical conflgurations jSl] and to the thermodynamics of fluids in porous 
media |35] . 

The test geometry we choose here is the case where the particle i3 is a single hard 
sphere S of radius i?s immersed in a pure hard-sphere fluid with radius R and density 
p. The change in grand potential Af2 due to the insertion of the sphere S is obtained 
by minimizing the density functional Eq. with either $wb or $wbii- From the 
equilibrium density proflle po(r) one can calculate Afi = r2[po(r)] — ^[p(r) = Pbuik]- 
By repeating the calculation for different values of i?s the function Af2(i?s) is obtained 
numerically. 

On the other hand, we have the morphometric prediction for Afi, Eq. ((201) • order 
to evaluate the morphometric solvation free energy it is most convenient to calculate 
the geometrical measures and the thermodynamic coefficients at the surface at which 
the density proflle po(r) jumps discontinuously to zero. This surface is parallel to the 
physical wall of S at normal distance R. Note that it is actually the parallel surface and 
not the physical surface that enters in the external potential Vcy±{r) in Eq. that S 
exerts on the fluid. In terms of the parallel surface, which is simply a sphere with radius 
Rs + i?, the morphometric form, Eq. ()2m) . reads 
4 

An(i?s) =p-n{Rs + Rf + (TA-K{R^ + Rf + KA-n{R^ + R) + RA-K. (21) 

The extraction of the thermodynamic coefficients p, a, /t and R from the values 
AVt{Rs) obtained by minimization of the density functional is therefore achieved by 
fltting Eq. (j2T|l to the numerical DFT data for different values of -Rs- This flt was 
performed for the data from $vvb and $wbii in the range R^ G [2i?, lOi?] for various 
values of the packing fractions of the fluid. Indeed, we flnd the assumption made by 
Eq. (PT|) on the i^g-dependence of Afi clearly conflrmed and in accordance with previous 
results jlHj. We shall come back to our results later on, referring to them as obtained 
via the "minimization route". 

In virtue of its applicability to mixtures, FMT also provides analytical expressions 
for the thermodynamic coefficients p, a, k and R. For their derivation, we follow the 
ideas of Refs. [221 EE] • In Ref. [221 the curvature dependence of the excess surface grand 
potential and contact density of a hard-sphere fluid in contact with hard curved walls 
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was studied. The basic idea is to consider a binary bulk mixture consisting of a hard- 
sphere fluid with radius R and packing fraction 77 and a single sphere 5, which is the 
second component at infinite dilution, i.e. ps 0. Af2 is then obtained as the excess 
chemical potential n'^. ^f^ can be calculated as the derivative of the mixture excess 
free energy density with respect to ps- Using any of the above FMT expressions for the 
excess free energy $ we find 

Pe^O 

Note that all vectorial contributions in $ vanish in the uniform bulk. Due to the 
limit Ps the partial derivatives of $ are evaluated for the solvent only, i.e. a one- 
component uniform fluid with radius R and packing fraction rj. 

A comparison of Eqs. ()21|) and (j^^ allows to identify the thermodynamic 
coefficients, calculated for the parallel surface, with certain linear combinations of the 
partial derivatives of $. We find that 

Pp 

(3a 











9$ 
dn2 




1 9$ 

Ait dni 




1 9$ 


R a$ 


47r driQ 


An drii 



pK = —- 2R— + R 



,2 



(23) 



dn 



3 



dn2 3 dns 

The relation for the pressure is precisely the scaled particle relation, Eq. ()10|). In 
the following, we refer to the above analytical results for the thermodynamic coefficients 
as the outcome of the "bulk route". 

We give the explicit results for the coefficients only for the case of the new excess 
free energy density $wbii: 

/3pwBii _ I + v + v'^ - r]^ 

p (1-^)^ 

pawBii _ 1 + 2r/ + 8r/^ - 5r/^ _ ln(l - r/) 

RP 3(1 3r/ ' 

pKwBU _ 4 - 10// + 20r/2 - 8r/3 4 ln(l - r]) ^ ' 

R^p 3(1 -r/)3 ^ 3^^ ' 

/?^WBii _ -4 + 1177-13772 + 4773 41n(l-7/) 

~ 3(1 -7/)3 3^^ ■ 

We emphasize that the pressure pweii is precisely the quasi-exact CS expression. 
This is not a trivial fact but rather a consequence of the construction of the novel 
mixture equation of state Eq. (|T5|) [Tlj . In contrast, the original White-Bear version of 
FMT, which is based on a different mixture generalization of the CS pressure 
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Figure 1. Results for the four thermodynamic coefficients p, a, k, and R of the hard- 
sphere fluid are shown as obtained from the new excess free energy density $\yBII' 
Eq. H16|l . The analytical expressions given in Eqs. I|24(l are denoted by the lines, while 
the results from the minimization route arc given by the symbols. 77 is the packing 
fraction. At w 0.494 the hard-sphere fluid freezes. 

does not possess this feature of self-consistency, i.e. the derivative of $wb with respect 
to ris does not yield exactly the CS pressure [TTj . 

We now compare the results from the bulk and minimization routes as obtained 
for the different versions $wb and $wbii of FMT. In Fig. Q we show our results for the 
thermodynamic coefficients calculated from the new functional $wbii- The agreement 
for p is perfect by construction of the equation of state, Eq. (fT3j). and very good for 
the surface tension a. Note that from a comparison with simulation data o"wbii was 
shown previously to be of high accuracy for intermediate and high packing fractions of 
the hard-sphere solvent. At low packing fractions, however, we found a small deviation 
from the exact low density limit of a Only for the bending rigidities k and R a 

slight inconsistency between the bulk and the minimization route appears. However, 
this inconsistency remains below 1% at high values of rj and we conjecture from the very 
good agreement of cwbii with simulation data that also kwbii and kwbii deliver accurate 
expressions for the hard-sphere fluid thermodynamic coefficients. With Eqs. (I24|l we 
have obtained a set of analytical expressions for the thermodynamic coefficients that 
are more accurate than previous suggestions, namely the results calculated from the 
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Figure 2. Various results for tlie four thermodynamic coeffieients, p, cr, k, and k. 
Shown are the differences of these coefficients obtained by various routes and theories 
to the analytical WBII results, Eqs. (|24|) . cf. the lines in Fig. ^ The symbols (except 
the crosses) denote the WBII results from the minimization route. For comparison, 
we also show results from the original White-Bear version: those from the bulk route 
are plotted as dashed lines while the crosses denote the outcome of the minimization 
route. 



original White-Bear version or those from Rosenfeld's DFT. 

As an illustration, we plot in Fig. |21 the difference of various results for the four 
thermodynamic coefficients, p, a, n, and K from the analytical expressions Eqs. ^I^i of 
the WBII version. Again, we find a high degree of self- consistency of the new functional 
$WBii (symbols, except the crosses in Fig. In contrast, the inconsistency of the 
original White-Bear version, which can be seen by the distance between the dashed line 
and the crosses in Fig. |21 is considerably larger and appears even for the pressure. The 
analytical expressions derived in the bulk route from $wb are therefore of a lower quality 
than Eqs. ()24|). which manifests itself also in their poorer agreement with simulations 
|14j . The good agreement between the results from the minimization route for the two 
versions of FMT is a direct consequence of the good agreement in the corresponding 
density profiles. This observation can be rationalized by noting that the contact value 
of the density profile at a planar wall coincides for both version of FMT as a result of 
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the same bulk equation of state. 

We do not include the results from Rosenfeld's functional $rf in Fig.|21as this would 
require to extend the range of the vertical axis considerably and therefore obscure the 
examination of consistency of the White-Bear versions. When calculated from $rf; 
however, the pressure follows the PY compressibility result which quantitatively differs 
from simulations so that the analytical expressions from the bulk route only yield a 
qualitative description of the thermodynamic coefficients. Surprisingly, the agreement 
between the bulk and the minimization route is comparable to that of $wb except 
for the pressure where $rf is consistent jHZj. Intuitively, one might expect a better 
agreement for than for $wbii because of the self-consistency on the level of the 
pressure. This interesting feature will be encountered again in the following where we 
consider the contact density at a curved hard wall. 

The contact density of the hard-sphere fluid at a hard wall is connected to the 
normal derivative of the grand potential VL j2H]- The case of interest here is again a 
hard-sphere fluid (radius i?, packing fraction rj) around a sphere S with radius Rs- 
We determine the grand potential from the density functional fi[po(r)]- The normal 
derivative of VL reduces due to the symmetry to a derivative with respect to Rs at 
constant chemical potential, which is then calculated as 

The flrst integral vanishes due to the equilibrium condition for /3o(r), i.e. SQ/6p = 0. 
The derivative of the external potential gives rise to a 6-peak at the location of the 
parallel wall, and one flnds pS] 
(90 

P^ = MRs + R?Pc (26) 

where pc is the contact value of the density of the fluid at the sphere S. Using the 
morphometric form Eq. (PT|) for AQ, the grand potential Q of the fluid containing the 
sphere S is fl{Rs) = —pVtot + AQ{Rs), where V^tot is the total volume of the system. In 
the thermodynamic limit Vtot — > C)0. If the morphometric form is inserted into Eq. (|26|l 
one obtains the contact density pc 

For Rs —* oo the planar wall contact theorem p^ = [3p is recovered and for flnite values 
of Rs the contact density is lowered. 

We show results for the contact density pc of a hard-sphere fluid with packing 
fraction rj = 0.4 as a function of the inverse of the radius Rs in Fig. IHl The 
symbols are the contact densities obtained from the density proflles which we calculated 
by minimizing numerically the functional Eq. (^. The lines in Fig. IHl show the 
morphometric prediction according to Eq. ^I7\i with the analytical expressions for the 
thermodynamic coefficients from the different versions of FMT (cf. the bulk route above). 
The flrst observation we make is that the numerical results from $wb and $wbii are 
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R/{R, + R) 

Figure 3. Contact value pc of the density of a hard-sphere fluid (radius R, packing 
fraction 77 = 0.4) at a single sphere with radius Rs- We show results obtained from the 
excess free energy densities $rf, *&wb and <I>wbii, respectively. We compare results 
from the numerical minimization of the density functional (symbols) , Eq. , with the 
morphomctric prediction (lines) according to Eq. (|77|l . 

nearly indistinguishable and indeed in the planar wall limit, Rs ^ 00, the data coincide 
by construction of the functionals. Taking this fact into account it is understandable 
that the results for pc for finite values of Rs are very similar. The numerical data for 
$RF tends towards the PY pressure for i?s — ^ 00 which is known to overestimate the 
actual pressure in the hard-sphere fluid for sufficiently high values of rj. In the limit 
Rs ^ (point-like object) the data from the three versions of FMT coincide. 

Comparing with the analytical prediction from morphometry, we find very good 
agreement between the results from $wbii over the whole interval of Rs. Only for very 
small radii Rs a slight deviation is visible. Therefore, the new functional improves upon 
the results of the original White Bear version of FMT, which performs well at small 
radii, but produces a small error at large values of Rs due to the the inconsistency of 
$WB in the pressure. 
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As mentioned above for the thermodynamic coefficients, a moderate agreement 
in the case of Rosenfeld's FMT is also observed for the contact density. While the 
approach is consistent for large values of Rs by construction of the functional, in the 
range of smaller values of Rs a deviation is clearly visible. This behavior has been 
observed previously jHS]. We find this fact remarkable because it shows that, from the 
point of view of self-consistency, the new hard-sphere mixture equation of state Eq. ()15p 
is better suited for an implementation within FMT than the PY mixture equation of 
state Eq. (|T3|) itself. This is even more surprising as the latter is characterized by full 
consistency for mixtures on the level of the pressure. 

5. Conclusion 

In this work, we have introduced a new density functional for hard-sphere mixtures 
which, in the spirit of the original White-Bear version, incorporates the quasi-exact CS 
equation of state within the framework of FMT. While the original White-Bear version 
is based on the well-known BMCSL equation of state, the White-Bear version Mark II, 
presented here, is derived from an improved mixture generalization of the CS equation 
of state recently introduced by us ^3] . The new functional WBII, besides having all the 
good properties of the original White-Bear version, improves the level of self-consistency. 
The level of consistency of the WBII version for the pure hard-sphere fluid is examined 
in the context of morphological thermodynamics. Our study reveals that, beside the 
improved consistency of the pressure, also the consistency of the surface tension a and 
the bending rigidities k, and k is clearly improved. Supported by a previous comparison 
to simulation data for the surface tension ^3] we can argue that the thermodynamic 
quantities derived from the new functional WBII are on a par with simulations. 

We have presented evidence that in the case of the pure hard-sphere fluid the degree 
of self-consistency of the WBII version is even higher than that of Rosenfeld's original 
FMT. This is a remarkable finding as Rosenfeld's FMT is by construction fully (i.e. 
for an arbitrary mixture) consistent on the level of the pressure. Apparently, this fact 
does not translate into a high level of consistency for other thermodynamic quantities 
(such as surface tension and bending rigidities) of the pure fluid. We conclude that, in 
what concerns the pure hard-sphere fluid, the recent mixture generalization of the CS 
equation of state Eq. (|T5|l is even better suited as a staring point for FMT than the PY 
compressibility mixture equation, underlying Rosenfeld's FMT. Note further that the 
PY mixture equation of state deviates significantly from simulations at sufficiently high 
densities. 

In conclusion, with the WBII version we have constructed a new hard-sphere 
functional based on the CS pressure which improves upon the original White-Bear 
functional. Although the differences between the density proffies in simple geometries 
resulting from minimization of the functionals are small, the increased self-consistency 
of the WBII version proves crucial for analytical calculations within the context of 
morphological thermodynamics. 
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Our considerations are based on thermodynamic arguments and result in a change 
of the dependency of the free energy density $ on the weighted density n^. There are 
several other developments in FMT that were mainly concerned with improving the 
performance of FMT in highly confined geometries. These studies suggest to change the 
dependency of $ on n2 and n2 |2(H l2T] or to introduce new tensorial weighted densities 
[221 • ■'■t is worth pointing out that these improvements concerning the description of 
hard-sphere fluids in highly confined geometries are straightforwardly combined with 
the improvements on thermodynamics presented here. 
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